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VO . ABSTRACT 

<N 

We study various classical solutions of the baby- Skyrmion model in (2 + 1) dimen- 
sions. We point out the existence of higher energy states interpret them as resonances 
of Skyrmions and anti-Skyrmions and study their decays. Most of the discussion in- 
volves a highly exited Skyrmion-like state with winding number one which decays into 
an ordinary Skyrmion and a Skyrmion-anti-Skyrmion pair. We also study wave-like 
solutions of the model and show that some of such solutions can be constructed from 
' the solutions of the sine-Gordon equation. We also show that the baby-Skyrmion has 

0^ ■ non-topological stationary solutions. We study their interactions with Skyrmions. 

In previous papers by two of us (BP and WJZ) [1-3] some hedgehog-like solutions 
of the so called baby-Skyrmion model were studied. It was shown there that the model 
has soliton-like topologically stable static solutions (called baby-Skyrmions) and that 
these solitons can form bound states. The interaction between the solitons was studied 

H 

in detail and it was shown that the long distance force between 2 baby-Skyrmions 
depends on their relative orientation. 

To construct these soliton solutions, one must use a radially symmetric ansatz 
(hedgehog configuration) and reduce the equation for the soliton to an ordinary dif- 
ferential equation. This equation admits more solutions than those described in [2] 
and [3], and, as we will show, they correspond to exited states or resonances made 
out of both Skyrmions and anti-Skyrmions. We will also show that the model admits 
some solutions in the form of non-linear waves. 

The (2 + l)-dimensional baby-Skyrmion field theory model is described by the La- 
grangian density 

L = F w (^djd a $- ^(djx d $)(d a $x 3p$) - /x 2 (l - n$)). (1.1) 
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Here <j>= (4>\,4>2,4>z) denotes a triplet of scalar real fields which satisfy the constraint 
(A 2 = 1; (d a d a = dtd 1 - d l d i ). As mentioned in [1-3] the first term in (1.1) is the 
familiar Lagrangian density of the pure S 2 a model. The second term, fourth order 
in derivatives, is the (2+1) dimensional analogue of the Skyrme-term of the three- 
dimensional Skyrme-model ' 4 l The last term is often referred to as a potential term. 
The last two terms in the Lagrangian (1.1) are added to guarantee the stability of a 
Skyrmion 

The vector n = (0,0,1) and hence the potential term violates the C>(3)-rotational 
iso-invariance of the theory. The state 3 = 1 is the vacuum state of the theory. 
As in [1,2] we fix our units of energy and length by setting F v = k = 1 and choose 
/i 2 = 0.1 for our numerical calculations. The choice n 2 = 0.1 sets the scale of the energy 
distribution for a basic Skyrmion. 

As usual we are interested mainly in field configurations for which the potential 
energy at infinity vanishes as only they can describe field configurations with finite 
total energy. Therefore we look for solutions of the equation of motion for the field $ 
which satisfies 



for all t. This condition formally compactifies the physical space to a 2-sphere S 2 h 
and so all maps from S 2 h to Sf so are characterised by the integer-valued degree of this 
map (the topological charge). 

The analytical formula for this degree is 



This degree is a homotopy invariant of the field <j> and so it is conserved during the 
time evolution. 

The Euler-Lagrange equation for the Lagrangian L (1.1) is 



One simple solution of (1.4) is given by ${x,t) = ft. This solution is of degree 
zero and describes the vacuum configuration. Another simple solution of (1.4) is 
evidently $ = -n. It is also of degree zero and may be considered as a false vacuum 
configuration. 

2. Static Skyrmions solutions 

Some static solutions of the equation of motion (1.3) were discussed in [1,2]. An 
important class of static solutions of the equation of motion consists of fields which 
are invariant under the group of simultaneous spatial rotations by an angle a e [0, 2n] 
and iso-rotations by -na, where n is a non-zero integer. Such fields are of the form 



lim^^oo^xj) = ft 



(1.2) 




(1.3) 



d a ($ x d a 4>- d^id^cf- $x d a $)) = n 2 (fxn. 



(1.4) 
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(sin f(r) cos{n9) \ 
sinf(r)sin(nO) , (2.1) 
cosf(r) J 

where (r, 0) are polar coordinates in the (x, y)-plane and / is a function which satisfies 
certain boundary conditions which will be specified below. Such fields are analogues 
of the hedgehog field of the Skyrme model and were studied in [3] for different values 
ofu 2 . 

The function /(r), the analogue of the profile function of the Skyrme model, has 
to satisfy 

/(0) = mi,meZ (2.2) 

for the field (2.1) to be regular at the origin. To satisfy the boundary condition (1.2) 
we set 



limr^cofir) = 0. (2.3) 

Then solutions of the equation of motion which satisfy (2.2) and (2.3) describe fields 
for which the total energy is finite. Moreover, the degree of the fields (2.1) is 

deg[<j>\ = [cos /(oo)- cos f '(())]- . (2.4) 

For fields which satisfy the boundary conditions (2.2) and (2.3) and which thus 
correspond to finite energy configurations, we get from (2.4) 

de ff M = [l-(-lH-. (2.5) 

The fields of the form (2.1) which are stationary points of the static energy func- 
tional V, the time independent part of L in (1.1), must satisfy the Euler-Lagrange 
equation for / 



n 2 sin 2 f , „,, . n 2 sin 2 f n 2 f sin f cos f . „, n 2 sin f cos f , „ 
(r + -)/" +1 5— + — /' - ru 2 sin] = 0. 2.6 

In [1,2] it was shown that the solutions of (2.6) for n = 1, 2 and m = 1, correspond to 
the absolute minima of the energy functional of degree 1 and 2 respectively. Any field 
obtained by translating and iso-rotating the solution corresponding to n = m = 1 was 
called a baby Skyrmion. As this field configuration has a topological charge deg[$] = 1 
we call it a "baryon" and denote it by the symbol B. A solution corresponding to 
n = -1 and m = 1 is then an anti-Skyrmion or an antibaryon B. 

Clearly, there exist also solutions of the equation of motion which satisfy the 
boundary conditions (2.2) at the origin but which, at infinity, behave as 

limr^oofir) = Itt. 

(2.7) 

I e z. 
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For such fields, if I is odd, the solutions (2.6) differ from B-Skyrmions and belong 
to a class of solutions with infinite energies. 

It is also worth mentioning that the asymptotic behaviour of these oo energy 
Skyrmion-like solutions differs from the asymptotic behaviour of £?-Skyrmions. Con- 
sider, e.g., the case m = 1,1 = 0,n = 1 (B-Skyrmion). Then the asymptotic behaviour 
of this solution is given by [2] : 



fW-K^r)-^^— e~» r (2.8) 

where K 1 (x) is the modified Bessel function. On the other hand the asymptotic 
behaviour of the oo energy anti-Skyrmion-like solution corresponding to (n = l,m = 
2, I = 1) is given by 

C 

/( r )~ 7r+ s inO + 7), (2.9) 

where C and 7 are constants. As the energy of these solutions is infinite we will not 
discuss them further in this paper and so we will always assume that I = in what 
follows. 




Figure l.a : Energy and Topological Figure l.b : Energy 

charge density for the hedgehog solu- charge density for tl 

tion : n=I, m=f tion : n=2, m=f 



We have solved (2.6) for different values of n and m using a shooting method 
with the appropriate boundary condition (2.2) (2.3) . We have determined the profile 
function f(r) and the energy and topological density profiles for each of these solutions. 
The total energies for our solutions are as follows: 
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1.5642 


5.011 
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2.9359 


7.725 


14.185 


22.233 
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4.4698 


10.555 


18.350 


27.819 


4 


6.1145 


13.465 


22.633 


33.395 



When m = 1, the topological charge is given by n, and each configuration is a 
superposition of n Skyrmions. We know from [2] that only the first two are stable. 
Figure 1, shows the profiles of the energy and of the topological density of the states 
n = 1 and n = 2. 



When to = 2, the topological charge is and, as can be seen from Figure 2, these 
configurations correspond to a superposition of n Skyrmions and n anti-Skyrmions 
where the Skyrmions form a ring surrounding the anti-Skyrmion at the centre. We 
have integrated separately both the positive and negative parts of the topological 
charge, and have found them to be n and -n respectively, thus justifying our inter- 
pretation. 




Figure 2. a : Energy and Topological Figure 2.b : Energy 

charge density for the hedgehog solu- charge density for tl 

tion : n=l, m=2 tion : n=2, m=2 



When to = 3, the topological charge is given by n, and each configuration is a 
superposition of 2n Skyrmions and n anti-Skyrmions (Fig. 3). The configuration is 
made of 3 layers, with n Skyrmions at the centre, n anti-Skyrmions in the middle, 
and n Skyrmions in the outside ring. 
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Figure 3. a : Energy and Topological 
charge density for the hedgehog solu- 
tion : n=l, m=3 



Figure 3.b : Energ 
charge density for 
tion : n=2, m=3 



Note that when m is larger than 1 the energy of the configuration is larger than the 
energy of nm baby Skyrmions. This indicates that such configurations are unstable 
and we will now analyse their decay modes. 

3. Exited meson Mi )2 (n = l;m = 2) 

Let us look first at the field which corresponds to n = 1, m = 2. From (2.4) we 
see that the topological charge of this field configuration is zero. So we can call this 
state, a "meson" or a coherent "meson cloud". We note from Figure 2. a that, like a B 
Skyrmion, this solution corresponds to a radially symmetric extended configuration 
with the maximum of the energy density at the origin (r = 0). 

As the total energy (5.011) exceeds the sum of the masses of a B-Skyrmion and a 
B-Skyrmion (2 x 1.5642), this configuration must be unstable. 

The distribution of the topological charge density shows more structure (Fig 2. a.) 
We note that the topological charge density is negative for small r and that it changes 
sign at r = r cr = 2. At this point the profile function /(r cr ) = n, thus the solution 
corresponding to the M ij2 state looks like a Skyrmion surrounded by an anti-Skyrmion 
field. Of course the total topological charge is zero. We see that the solution still looks 
like an extended but localised configuration. However, although the topological charge 
density suggests that the Skyrmion is at the origin the profile function there is given by 
/ = 2tt and so, from this point of view, resembles more the vacuum than a Skyrmion. 
This suggests a possible interpretation of the Mi, 2 state in terms of Skyrmions and 
anti- Skyrmions of the usual or oo energy type. Whatever the interpretation, the state 
is mesonic in nature. Moreover, looking at the field configuration we note that most 
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of its changes takes place around those points in the (x,y) plane where r = r cr i.e. 
where / = n. It is remarkable that everywhere along this circle 4> 3 = -1. This circle 
is actually the region of instability when the solution M lt2 is excited by a non-radial 
perturbation of small amplitude. 

In fact, it is not difficult to demonstrate that at each point on the circle r = r cr 
one can create "hedgehog-like" extended objects using only infinitesimal perturba- 
tions. Of course, if these extended objects were to correspond to Skyrmions and 
anti-Skyrmions then to conserve the topological charge they will have to be created 
in pairs and we would expect them to appear as soon as we perturb the initial con- 
figuration. To check for such a behaviour we have perturbed the initial configuration 
corresponding to the Mi^-state. Our perturbation was in the form of small excess 
of kinetic energy centered around two points in the (x,y)-plane chosen symmetrically 
with respect to r = 0. We have found that, indeed, the state split into a Skyrmion and 
anti-Skyrmion pair of the B-type. Their relative orientation was such that the force 
between them was repulsive. After their creation the Skyrmion and the anti-Skyrmion 
moved in opposite directions from the centre. 

We then experimented with perturbing the initial state by different perturbations 
and we observed different decay modes. When the applied perturbation was not 
symmetric with respect to r = (but was close to being symmetrical) the M lj2 state 
decayed into a Skyrmion and an anti-Skyrmion, which then rotated in their internal 
space so that their relative orientation made them to attract each other. They then 
collided into each other and decayed into waves. 

We thus conclude that the solution Mi >2 of (2.6) is indeed a saddle point in the 
space of field configurations. The state can be thought of as a resonance of a Skyrmion 
and an anti-Skyrmion and it has different decay modes when perturbed. It decays 
either into a BB-p&ir or into light "mesons" . 

It would be interesting to see whether it is possible to create the M 12 state in a 
head-on collision of a Skyrmion and an anti-Skyrmion. We are planning to come back 
to this question in a future paper. 

4. Exited Bi <3 baryon (n = l,m = 3) 

This is a baryon-type state as its topological charge is one. The profiles of the 
energy and of the topological charge densities for this state are shown in Figure 3a. 
Looking at the density of the topological charge we see that this state may be consid- 
ered as the usual B-Skyrmion surrounded by an anti-Skyrmion and a further Skyrmion 
ring. The configuration has an energy larger than the energy of its constituents and 
is thus unstable. 

We have performed several simulations looking at the decay products of this state. 
When we used a perturbation symmetric with respect to the origin the configuration 
decayed in 2 Skyrmions and 1 anti-Skyrmion, the anti-Skyrmion staying at the ori- 
gin, while the two Skyrmions moved in opposite directions. When the perturbation 
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was not symmetric, the Skyrmions and the anti-Skyrmion were able to change their 
relative orientations and one of the Skyrmions collided with the anti-Skyrmion and 
decayed into waves. 



5. Mesons made from dibaryons and anti-dibaryons 

Let us now discuss various properties and decay modes of exited meson-like states 
made out of dibaryons (a bound state of 2 B-Skyrmions) and antidibaryons. We will 
concentrate our attention on the n = 2 and m = 2 state but our discussion generalises 
easily to other states. The topological charge of the (n = 2, m = 2) configuration (ie 
its baryon number) is zero, so this configuration is a meson-like state from the point 
of view of our classification. As n = 2 the configuration looks like a dibaryon near the 
origin, i.e. near r = 0. From Figure 2.b, we see that it corresponds to a dibaryon at the 
origin surrounded by an anti-dibaryon ring. The border between these two regions of 
opposite topological charge is, again, very well defined and is situated along the circle 
of radius r cr = 3, (f(r cr ) =ir). So this ring is again the region of instability. Applying 
different types of perturbations to our (n = 2, m = 2) solution, we have observed three 
different decay modes for this state: 



The pictures of the energy density for the first decay modes is shown in Figure 4. 



Figure 4 : Topological charge density 
for the M 2 , 2 -» B+B+B+B decay mode. 

6. Other exited mesonic and baryonic states 

i) Let us look first at the state (n = l,m = 4). This state, again, is mesonic with 
deg[$] = 0. The energy and topological charge profiles are shown in Figure 5. 




waves 



(5.1) 
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The state is more complicated as its energy density exhibits additional maxima 
and minima. In fact, the state looks as if it were a coherent state of Skyrmions 
and anti-Skyrmions, with rings of different radia occupied by fields of alternating 
topological charge. This is clearly seen from the topological charge density plots; 
moreover, the total topological charge in each ring is +1,-1,+1 and -1 respectively 
(going out from r = 0). Thus we denote this state as M lj4 meson. 




Figure 5. a : Energy and Topological Figure 5.b : Energy 

charge density for the hedgehog solu- charge density for tl 

tion : n=l, m=4 tion : n=2, m=4 




Figure 6 : Topological charge density 
for the M 2 , 3 -> B+B+B+B+B+B decay 
mode. 
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ii) Another interesting state is that of M 2j3 (n = 2,m = 3). Its topological charge is 
2. Hence this state can be thought of as an exited state of the dibaryon M 2 ,i. 
Its energy is clearly quite large and the state represents a coherent mixture of 4 
Skyrmions and 2 anti-Syrmions and so it can decay into different channels. One 
example of its decay mode is shown in Figure 6. We see in this picture that the 
decay products consist of two outgoing Skyrmions and two anti-Skyrmions with 
the remaining 2 Skyrmions at rest at the origin. The fact that the decay products 
involve one di-Skyrmion at the origin (slightly excited) is due to the attraction 
between two Skyrmions, as mass(M 2j i) < 2mass(i?), see our table in section 2. 

iii) Another interesting state is M 2j4 which corresponds to (n = 2, m = 4) and, as such, is 
a highly excited meson state. Its energy and topological charge profiles are shown 
in Figure 5b. 

Clearly, the list of new solutions may be continued further by taking larger values 
of n and m. Of course, having calculated some of their profile functions we see that as 
one increases m their energies increase (quite rapidly). All these higher energy states 
can be treated as coherent states of some number of Skyrmions and anti-Skyrmions 
and are unstable. Under suitable perturbations they will decay into various channels 
involving Skyrmions and anti-Skyrmions with some Skyrmions and anti-Skyrmions 
annihilating into pure waves. 

7. Plane wave solutions 

Let us now look at wave-like solutions of our model. Recall that such wave- 
like solutions have already been studied in ^ for a slightly different version 
of a "skyrme-like" (2+1) dimensional model (the potential term of that model was 
different.) 

What type of waves can we find for the Lagrangian in the form (1.1) ? To look for 
plane wave-like solutions, we seek solutions which do not depend on one variable, say 
y. However, as soon as we impose this condition we note that the Skyrme-like term 
vanishes for these field configurations. So the discussion is not that difficult and can 
be performed analytically. 

First, we look for solutions of the equation of motion for the field $ in the form: 

4> = (sinfcosip, sinfsinip, cos /), (7.1) 
where / = f(x,t) and V = ip(x,t). In terms of / and ip the Lagrangian takes the form 

L = (1/2) J d 2 x (djd^f + ((1 - cos(2f))/2)d^d^ - ^(1 - cos(f)). (7.2) 

To go further we consider the case when the phase ip is constant. In this case the 
field / satisfies 

d^f + i? sin(f) = 0, (7.3) 
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which is the sine-Gordon equation. So we see that when tp = const the wave solutions 
of the Lagrangian (1.1) are given by the solutions of the sine-Gordon equation. 

Moreover, solutions of (7.3) with small amplitude may be considered as ordinary 
plane waves with the dispersion relation given by 



u z = n z + k\ (7.4) 

Of all finite energy solutions of (7.3) (one-dimensional case) the most important, 
and perhaps the best studied, is the solution of the kink type: 

f(x, t) = 4 atan(exp [ - p( £z£lZJ^ )] ) . (7.5) 

When translated to our case we note that when x changes from -00 to +00, the 
vector $ moves along the meridional cross-section of the Sf so -sphere and returns to 
the same point. Moreover, this is true for any fixed time t. 

Another solution of (7.3) , called breather, is also well known and has been studied 
by many people. Its form is ' 8 ' 

/(,,t) = 4atan( ( 1 -^ 1/2 J^'lf J (7.6) 
V lo cosh{{\ — uj z ) L ' z (x — .t )) / 

where u is the frequency of the internal breather oscillations. 

These 2 types of solitonic waves are infinite front lines (the solutions depend only 
on one variable, x). Unfortunately, once imbedded into S 2 , the target space for our 
model, the extra degrees of freedom make these waves unstable. For the kink solution 
this is not surprising as the "loop" around the meridian can easily "slip" on one side 
of the sphere, thus decreasing the potential energy. 

To see this, let us take the configuration 

(cosa sin/ \ 
sinacosa (1 - cos/) , (7.7) 
1 - cos 2 a (1 - cos/) / 

where / is given by (7.5) with v set to 0, and where a is a function which depends 
only on y and which goes to when y goes to infinity. This configuration describes 
the kink (7.5) (7.1) perturbed locally in x and y, and it "displaces" the kink from the 
meridian of the 2-sphere onto a loop of a smaller radius on S 2 . The energy density 
for this static configuration is given by 

E = [ dxdy - |cos 2 a;/ 2 + a 2 (sin 2 asin 2 f + (1 - cos/) 2 ) 

J 12 (7.8) 
+fc 2 / 2 a 2 (l-COS/) 2 COS 2 a] + /i 2 COS 2 a(l - COS/) . 
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To prove that for some appropriate choice of a, the energy decreases, we compute the 
change of energy induced by a non-zero a 



SE =E(a = 0) - E(a) 

= J dxdy 2sin 2 ^ [fi 2 (sm 2 a 

+ ^i 2 sin 2 a(l - cos/) 

> I dxdy 2sin 2 ^[M 2 (sin 2 a 




fc 2 a 2 COSa) - a 2 y (l + sin 2 a)] + ^ 2 sin 2 a(l - COS/) 



(7.9) 



where we have used the fact that for (7.5) f x = 2/ism(//2). Given a(y) satisfying the 
asymptotic behaviour imposed, we can always stretch it by performing the dilation 
y — > ay to make a y small enough to make the two negative terms in (7.9) smaller than 
the positive terms. This proves the instability of the sin-Gordon kink wave when 
imbedded into the S 2 model. 

We have performed some numerical simulations for both types of waves and have 
indeed observed their instability. In both cases, as soon as some region of the wave 
is perturbed, the wave collapses around this point, emitting radiation. The collapse 
front then propagates rapidly along the solitonic wave destroying it completely. 

In our previous work ^, we have studied the scattering properties of plane waves 
and skyrmions. The skyrmions studied in [7] were different (the potential was dif- 
ferent) so we have repeated our analysis for the model studied here and have found 
no major difference. Like in the previous case, the skyrmion absorbs a section of the 
wave and starts moving after the collision. 

We have also studied the scattering between a skyrmion and the sine-Gordon front 
waves. We observed that when the skyrmion is close to a sine-Gordon wave (breather 
or kink) it triggers the wave collapse and, as a result, there is no real scattering 
between these two objects. 

While performing simulations with large amplitude breather waves we have ob- 
served the formation of a radially symmetric breather-like soliton. This solution 
looks very similar to the pulsons observed in [9]. We discuss some of its properties 
in the next section. 

8. Non-topological Solitons 

In the previous section we have described our studies of waves in the baby-skyrmion 
model in which we observed that the plane wave solutions which are given by the 
solutions of the (1+1) dimensional sin-Gordon equation are unstable. 

However, we when we looked at the decay of some breather front waves we encoun- 
tered something rather unexpected. Instead of decaying into waves that dissipate like 
kink solutions they produced a radially symmetric breather-like field configuration 
which appeared to be relatively stable. By looking at the time evolution of this field 
configuration, as produced by our simulation, we were able to conclude that: 
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- the field configuration is radially symmetric. 

- up to a global rotation of S 2 , the solution "lives" in the <pi,<f>3 plane of the target 
space (S 2 ). 

To study this field configuration further we make the following ansatz: 

= (sm/(r,t),O,cos/(r,t)). (8.1) 
The lagrangian density then becomes 

L = ttJ drr(d t fd t f - d r fd r f - y 2 (l - COS(/)). (8.2) 

and the equation reduces to the radial sin-Gordon equation: 

ftt ~ frr - y + ^Sm(f ) = 0. (8.3) 

This equation has already been studied in [9] , where it was shown that it has time 
dependent solutions similar to a breather, but which radiate their energy and slowly 
die out. The authors of [9] decided to call such configurations pulsons. 

In [10] we have also shown that there exist stable time dependent solutions. The 
radial field configurations of (8.3) radiate relatively quickly when their amplitude 
of oscillation is relatively small, that is when the value of / never becomes larger 
than 7r/2 at the origin (these are the pulsons studied in [9].) When the amplitude of 
oscillation is larger than n/2 the breather-like configuration slowly radiates its energy 
and asymptotically reduces the amplitude of oscillation to n/2 and and settles at a 
period of oscillation T — 20.5 (when = O.l). By trial and error we have found that 



/(r,0) = 4atan(cexp(--^atan(^) 



(8.4) 



with K = 10 and C = tan(7r/8) is a good initial condition for this metastable pseudo- 
breather solution. 

As in the case of plane wave solutions we have to check that once they are imbedded 
into S 2 they are still stable. We have indeed checked this numerically and we have 
found that the solution described by (8.1) and (8.3) is indeed stable in the S 2 model. 

The energy of this new breather-like solution is given by 

E PB ~ 3.97 

which means that it is 2.5 heavier than the baby-skyrmion. Moreover, its topological 
charge density is identically zero but it has enough energy to decay into a skyrmion 
anti-skyrmion pair. In what follows we shall refer to these field configurations as 
pseudo-breathers. 



13 



In practice, it is very difficult to have a field configuration of a pseudo-breather. 
However, we can find approximate field configurations which still radiate energy and 
asymptotically reach the stable (or perhaps only metastable) configuration config- 
uration of the pseudo-breather. The excess of energy over the final configuration 
can then be seen as an excitation energy which is slowly radiated away. During any 
scattering process solitons tend to exchange or radiate some energy. When the ex- 
citation energy is large enough, the outgoing pseudo-breather-like configuration may 
have enough energy to evolve into the stable pseudo-breather field; otherwise, it ends 
up with less energy than the metastable configuration and progressively dies out. 

The scattering properties of pseudo-breathers are quite interesting. When the 
pseudo-breathers are imbedded into the baby-skyrmion model the field configurations 
have an extra degree of freedom corresponding to their orientation inside the 4>i,<h 
plane. When two pseudo-breathers are set at rest near each other, the force between 
them depends very much on their relative orientation: when they are parallel to each 
other and oscillate in phase, they attract each other, overlap and form a new structure 
which appears to be an exited pseudo-breather. This pseudo-breather then slowly 
radiates away its energy. The non-topological nature of pseudo-breathers means that 
they can indeed merge to form a new structure of the same type. 

If the two pseudo-breathers are anti parallel, i.e. if they oscillate completely out 
of phase, then the force between them is repulsive. When the two pseudo-breathers 
have a different orientation they slowly rotate themselves until they become parallel; 
then they move towards each other and form an exited pseudo-breather structure. 
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Figure 7. a : Pseudo-Breather scatter- Figure 7.b : Pseudo 

ing. Impact Parameter: 10, v = 0.2 ing. Impact Paramet 

When two pseudo-breathers are sent towards each other with some kinetic energy, 
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the scattering is more interesting. Depending on the initial speed or the scattering 
impact parameter, they either merge into a single pseudo-breather or they undergo 
a forward scattering. The details of these scattering properties are given in [10]. In 
Figure 7 we show two such trajectories. For both scatterings the initial speed was 
0.2, the only difference being the values of the impact parameter. 

9. Pseudo-Breather-Skyrmion Scattering 

As we have seen, the baby-skyrmion model has two different types of extended 
solutions. The skyrmions are topological solitons which are very stable, while the 
pseudo-breathers are time-periodic solitons which can slowly decay if they are per- 
turbed too much. It is very unusual to have a model that exhibits two such different 
stationary structures and so it is interesting to analyse how they interact with each 
other. 

Total Energy : t = 0.000000 



Min = 0.000009 Max = 0.165958 




Figure 8 : Skyrmion and pseudo-breather at rest. 

When a skyrmion and a pseudo-breather are put at rest next to each other the 
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overall interaction between them makes the skyrmion slowly move away from the 
pseudo-breather-like configuration while the pseudo-breather looses some of its energy 
faster than when it is placed there by itself. 

To scatter a skyrmion with a pseudo-breather we have placed the pseudo-breather 
soliton at rest, and we have sent the skyrmion towards it. We have performed this 
scattering for different orientations of the pseudo-breather, for different values of the 
impact parameter and for a few different speeds. In each case, the pseudo-breather 
was initially located at the origin while the skyrmion always started from (x ,y ) 
where xo is the initial position along the x axis and y is the impact parameter. The 
pseudo-breather was oscillating in the (<pi,<p2) plane along the direction (cos(a),sin(a)). 

In Figure 8 we show the initial condition corresponding to a baby-skyrmion next to 
a pseudo-breather. We note that the skyrmion is much more spiky than the breather. 

Our numerical results are summarised in the following 4 tables. 
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Table l.a : Impact parameter and speed dependance of the scattering angle (a = 0, x = -20) 
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Table l.b : Impact parameter and speed dependance of the scattering angle (a = w/2, x Q = -20 
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Table l.c : Impact parameter and speed dependance of the scattering angle (a = 7r/4, x n = -20 



v \ xO 


15.2 


16.2 


17.2 


18.2 


19.2 


20.2 


0.2 


38 


33 


33 


38 


31 


21 


0.4 


15 


16 


17 


20 


22 


19 



Table 2 : Scattering angle as a function of the initial distance (a = w/2, y = 2.5). 

The amount of energy lost by the pseudo-breather during the scattering is larger 
when the overlap between the skyrmion and the breather, both in time and space, 
increases. In some cases the pseudo-breather is completely destroyed by the collision. 
The oscillation of the pseudo-breather makes the interaction time dependent, and, as 
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a result, it is difficult to extract a simple pattern from the tables of scattering angles. 
10. Conclusions 

We have shown that the baby-skyrmion model has many interesting classical so- 
lutions in addition to the skyrmion solitons. The first class of solutions describe 
exited states of skyrmions and anti-skyrmions which are unstable with respect to 
perturbations. 

The second class of solutions involves non-topological stationary stationary field 
configurations which are periodic in time. They are relatively stable but can be 
destroyed if they are perturbed too much. 
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ABSTRACT 

We study a localised metastable breather-like soliton of the baby-Skyrmion model 
in (2 + 1) dimensions. We show the relative stability of that solution and study how 
it scatters with the skyrmion solitons. 

1. Introduction. 

In a previous papers , we have analysed the existence and some properties of waves 
in the so called baby-Skyrmion model. In that paper we announced that the model 
possesses a localised breather-like solitons. We are now studying this new solutions in 
more details, including its scattering properties with the other solitons of the model 
(the skyrmions). 

The (2 + l)-dimensional baby-Skyrmion field theory model [2-4] is described by the 
Lagrangian density 

L = F 7T (±dJd a $- ^(djxd $)(d a $xd p $)-fi 2 (l-n$)). (1.1) 

where $= (<j>i,<j>2,<j>3) denotes a triplet of scalar real fields satisfying the constraints 
$ 2 = 1 and (d a d a = d t d l - did 1 ). The first term in (1.1) is [2-4] the familiar Lagrangian 
density of the pure S 2 a model. The second term, fourth order in derivatives, is the 
(2+1) dimensional analogue of the Skyrme-term of the three-dimensional Skyrme- 
model ' 5 l The last term is a potential term necessary ^ to have stable topological 
solitons (skyrmion). 

The vector ft = (0,0,1) and hence the potential term violates the 0(3)-rotational 
iso-invariance of the theory. The state <f> z = 1 is the vacuum state of the theory. 
As in [4,3] we fix our units of energy and length by setting F n = k = 1 and choose 
H 2 = 0.1 for our numerical calculations. The choice n 2 = 0.1 sets the scale of the energy 
distribution for a basic Skyrmion. 
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The Euler-Lagrange equation for the Lagrangian L (1.1) is 

d a ($xd a $- d p $(d p $- $x d a $)) = n 2 $xn. (1.2) 

As shown in [3] and [2], the model admits topological solitons which we call baby- 
skyrmions. The topological charge of any field configuration is given by 

Q = -L J $■ (d^x d 2 $)d 2 x. (1.3) 

The simplest solution (Q=l) is obtained from the so called hedgehog ansatz 

sin f(r) cos(nQ) 

cf)(x) = I sin f (r) sin(n6) I , (1.4) 



cos 



f(r) 



where (r,6) are polar coordinates in the (z,t/)-plane and / is a function which satisfies 
the boundary conditions 

f(0) = mTr,m£Z (1.5) 

and 

linir^oofir) = 0. (1.6) 

The fields of the form (1.4) which are stationary points of the static energy functional, 
the time independent part of L in (1.1), must satisfy the Euler-Lagrange equation for 
/ 

n 2 sin 2 / „ n 2 sin 2 f n 2 /' sin f cos f , n 2 sin f cos / 2 . 
r + / + 1 + / rnsmf = 0. 1.7) 

The baby-skyrmions are stable solutions, and it was shown in [3] that they attract 
each other, and can form stable bounds states. 

2. Breather like solutions 

In [1] we have studied the existence of waves in the baby-skyrmion model. We 
have shown that there are plane wave solutions which are given by the solutions of 
the (1+1) dimensional sin-Gordon equation. 

Indeed, if we look for solutions of the equation of motion for the field ^in the form: 

4> = (sin f cos %p, sin f sinijj, cos /), (2-1) 
where / = f(x,t) and ip any constant, the Lagrangian takes the form 

L = (1/2) J d 2 x (djd^f -ii 2 {l- cos(f)). (2.2) 

and the field / satisfies 

a (i a"/ + /t 2 sm(/)=0. (2.3) 
2 



These one dimensional waves correspond to front waves when embeded into the 
(2+1) dimensional baby-skyrmion model. 

Unfortunately, these front waves solutions are unstable' 1 ^. Both the sin-Gordon 
kink, and the breather solutions decay as soon as they are perturbed. 

The suprise came from the decay of the breather front waves. Instead of decaying 
into waves that dissipates like the kink solutions, it produced a radially symmetric 
breather like solution which appeared to be relatively stable. By looking at the field 
configuration, produced by our simulation, we were able to draw two facts: 

- the solution is radialy symmetric. 

- up to a global rotation of S 2 , the solution "lives" in the <f>i,4>3 plane of the target 
space (S 2 ). 

We can thus make the following ansatz for the localised breather: 

$=(sinf(r,t),0,cosf(r,t)). (2.4) 
The largrangian density becomes 

L = it J drr(d t fd t f - d r fd r f -ix 2 (I- cos(f)). (2.5) 
and the equation reduces to the radial sin-Gordon equation: 

ftt -f rr - — +» 2 sin(f) = 0, (2.6) 
r 

This equation has already been studied in ' 7 ^, where is was shown that it has time 
dependant solutions simillar to a breather, but which radiates its energy and slowly 
dies out. They decided to call this solution a pulson. 

One of our first observation about this solution is that is radiates relatively quickly 
when its amplitude of ascillation is relatively small, that is when the value of / never 
becomes larger than tt/2. When the amplitude of ascilation is larger than 7r/2, then the 
breather progressively radiates its energy and assymptoticaly reduces the amplitude 
of oscillation to ir/2 and a period of oscilation T = 23 (when /t 2 = O.l). 

This suggest that there is a metastable solution of amplitude n/2. We have per- 
formed very long simulation with up to 20000 oscillation, and have observed that 
the amount of energy radiated by the pseudo-breather was decreasing steadily. The 
energy of the configuration assymptoticaly settled to E = 4.4. 

By trial and error, we have found that the extression 

f(r,0) = Cexp(-(r/K) 2 ) 
with K = 12 and C = 2.25, settles rapidly to the metastable solution. 

3. Breather-Skyrmion Scattering 
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Impact parameter and speed dependance of the scattering angle, a = 
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Impact parameter and speed dependance of the scattering angle, a = n/2 
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Impact parameter and speed dependance of the scattering angle, a = n/4 
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Scattering angle as a function of the initial distance, a = n/2,y0 = 2.5 



4. Breather-Breather Scattering 
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